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Abstract: 


In this paper we consider the effects of the isotopic composition of the primary galactic 
cosmic rays (GCR), nuclear fragmentation cross-sections, and isotopic-grid on the 
solution to transport models used for shielding studies. Satellite measurements are used to 
describe the isotopic composition of the GCR. For the nuclear interaction data-base and 
transport solution, we use the quantum multiple-scattering theory of nuclear 
fragmentation (QMSFRG) and high-charge and energy (HZETRN) transport code, 
respectively. The QMSFRG model is shown to accurately describe existing 
fragmentation data including proper description of the odd-even effects as function of the 
iso-spin dependence on the projectile nucleus. The principle finding of this study is that 
large errors (+100%) will occur in the mass-fluence spectra when comparing transport 
models that use a complete isotopic-grid (-170 ions) to ones that use a reduced isotopic- 
grid, for example the 59 ion-grid used in the HZETRN code in the past, however less 
significant errors (<+20%) occur in the elemental-fluence spectra. Because a complete 
isotopic-grid is readily handled on small computer workstations and is needed for several 
applications studying GCR propagation and scattering, it is reco mm ended that they be 
used for future GCR studies. 
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INTRODUCTION 

An important goal for NASA’s Space Radiation Health Program is to develop a 
predictive capability to predict the GCR fluence spectra to within a +25% accuracy 
(Annon. 1998). NASA has developed the HZETRN (high-charge and energy transport) 
code (Wilson, 1977; Wilson and Badavi, 1986; Wilson et al., 1991) as a science 
application and engineering design tool (Wilson et al., 1993) to be used in space radiation 
shielding studies. HZETRN has been validated in its ability to predict total dose and dose 
equivalent behind several materials in space to within +20% on multiple space missions 
in Earth orbit (Cucinotta et al., 2000a, Badhwar and Cucinotta, 2000, Badhwar et al., 
2001). However, interest in fluence-based approaches to risk assessment (Cucinotta et al., 
1995a, 1996a) suggests that more stringent tests of transport code accuracy be made, and 
the quantities dose and dose equivalent are deemed as necessary, but not sufficient tests 
of their accuracy. In this regard, we note that dose and dose equivalent are integral 
quantities that receive contributions from many GCR charge groups. There currently exist 
large uncertainties in biological response models for GCR (Cucinotta, et al., 2001) and 
methodologies to estimate health risks such that dose and dose equivalent may be 
insufficient as tests of transport code accuracy. The use of ion fluence as a basis for tests 
for accuracy provides for sufficient generality to ensure accuracy in GCR transport 
models including under the circumstances of revision of radiation quality factors or 
integration of alternative risk assessment approaches in the future. 

In the description of the transport of the galactic cosmic rays (GCR) in shielding 
materials or tissue, a common approximation is to consider only the elemental 
composition of the primary GCR and a reduced isotopic-grid for the secondary nuclei 
produced in nuclear fragmentation. In this paper we analyze the role of the isotopic 
dependence of the GCR primary composition and nuclear fragmentation in predicting the 
fluence of the GCR behind arbitrary shielding configurations. Our study is an important 
milestone in achieving NASA’s goal of accurate GCR transport codes, since for the first- 
time a complete isotopic-grid has been achieved in a GCR transport model and we 
document the error inherent in former approaches. Also, for applications that will 
consider radioactive isotopes produced in the atmosphere or shielding, our study provides 
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a useful tool to perform such analyses. Other applications where non-stable nuclei are 
considered are studies of the origin or the GCR where so-called cosmic-ray “clocks” 
consider the primary or secondary GCR with life-times on the order of confinement time 
in the galaxy (~1 M-yrs) (Yanasak et al., 1999). Several GCR “clock” nuclei including 
10 Be and 26 Al were not considered in the grid used in HZETRN in the past. Finally new 
data on the GCR near Mars are being collected by the MARIE experiment on the 
Odyssey spacecraft (Zeitlin et al., 2004), and our study provides an opportunity to begin 
new investigations on the accuracy of computational models used to describe the GCR. 

Historically the HZETRN code grew from a 29-ion isotopic grid used in the 1980’s and 
early 1990’s (Wilson, et al., 1991) to an extension to a 32-ion isotopic grid made in 1993 
in order to include all light ions (Cucinotta, 1993). Because of the limitations of random 
access memory present on the computer workstations of the early 1990’s, sensitivity 
studies were made for mono-energetic ion beams to study the minimum number of 
isotopes for convergence resulting in the use of 59-isotopic grid (Kim, et al., 1994), and 
all GCR studies since 1994 have used the 59-isotopic grid (Shinn, et al., 1994). However, 
there are several reasons to re-consider the use of the full isotopic-grid for GCR transport 
problems. First, the isotopic dependence of the primary GCR has not been considered in 
past shielding studies and may lead to errors in the description of both primary ion 
attenuation and secondary particle production including the role of high-energy neutron 
production from the many neutron rich species that occur. Secondly, the studies of Kim et 
al. (1994) used the NUCFRG2 model of fragmentation (Wilson et al., 1994), which does 
not provide a correct description of the even-odd effect observed in fragment production 
or of the projectile iso-spin dependence observed experimentally (Knott et al., 1996, 
1997; Zeitlin et al., 2001). Thirdly, the sensitivity studies made by Kim et al. (1994) used 
a “calibration” of the isotopic-grid to 56 Fe beams; however a larger isotopic grid occurs 
when all GCR projectile nuclei are considered. Fourth, the error in the range-energy and 
stopping powers that results from the use of a reduced isotopic grid, although expected to 
be small for large mass number, A»l, is an unnecessary one for transport calculations. 
Finally, the improved computational speed and memory available on current small 
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computer workstations, makes the inclusion of a complete isotopic-grid in the HZETRN 
code to be readily implemented at this time. 

In this paper the implementation of the HZETRN code to include the full isotopic 
dependence of the primary GCR is described. The physics of isotopic effects in GCR 
transport are described and the fragmentation parameters are a key component of this 
description. The quantum multiple scattering theory of nuclear fragmentation (QMSFRG) 
is used as the generator for fragmentation cross-sections used in our study. An empirical 
model of the isotopic composition of the primary GCR including its solar modulation is 
also described. For GCR problems an isotopic grid of 170 ions is identified and 
comparisons made to previous HZETRN results using the reduced-grid of 59 isotopes are 
made. The present code includes all of the abundant nuclei in the GCR environment with 
fluxes greater than about 10 /cm /yr and nuclei produced in fragmentation events with 
production cross-sections greater than about 1 mb. Several nuclei with smaller primary 
abundances or production cross-sections, which are of interest for scientific reasons are 
also included in the expanded HZETRN model described herein. The resultant code 
includes many neutron rich nuclei that have been ignored in the past with iso-spin 
components ranging from T-= +3/2 to T z =-3. 

ISOTOPIC COMPOSITION OF THE GCR 

NASA currently uses the GCR model of Badhwar and O’Neill (1992) to describe the 
elemental composition and energy spectra of the GCR including their modulation by the 
sun’s magnetic field. In this GCR representation, only the most abundant GCR nuclei is 
considered for each element and other isotopes of identical charge are counted as the 
abundant isotope. However, theoretical models and satellite measurements of the GCR 
have long considered the isotopic composition of the GCR and their modification through 
transport in interstellar space including estimating the primary nuclear composition at 
stellar sources (Parker, 1965, Webber et al., 1990a, Fields et al., 1994). The approach 
used here is to estimate an energy-independent isotopic fraction, fj from satellite 
measurements, which are constrained to obey the sum-rule 
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(1) <t>(Z,E) = 


where the left-hand side of eq. (1) is the elemental spectra from the Badhwar and O’Neill 
model and V ,/ y = 1 . Equation (1) is used herein as an initial estimate of the influence 

of the primary isotopic composition on GCR shielding calculations. Experimental studies 
have included measurements on the Pioneer, Voyager, and Ulysses spacecraft. A survey 
of such data (Hesse et al., 1991; Lukasiak, et al., 1993, 1995; Webber et al., 1985, 1990; 
Wiedenback et al, 1981, 1985) was made with the results shown in Table-1. 

The GCR path-length distribution represents the mean amount of inter-stellar and inter- 
planetary material intersected by cosmic rays prior to their arrival in the near-Earth 
environment. This distribution is currently estimated to vary between 3 and 20 g/cm with 
the inter-planetary gas approximately 90% H and 10% He. Variability in estimates 
between various studies arise because of differences in model nuclear fragmentation 
parameters, cosmological interaction terms, and factors related to the near-Earth or deep 
space measurements including energy and mass resolution of the detectors as well as 
position in the solar cycle. In most current models, the GCR path-length distribution is 
estimated to be velocity and rigidity dependent (i.e., energy dependent). Of note is that 
based on the solar modulation theory of Parker (1965), the inter-planetary portion of this 
distribution would be modulated over the solar cycle suggesting the form 

(2) /L(d>,R) = 2 

Inter -stellar ^ Inter -planetary (d>,R) 

where & is the solar modulation parameter and R (in A.U.) is the radial distance from the 
sun. It then follows that the amount of fragmentation in the inter-planetary media and the 
isotopic abundances are dependent on the modulation parameter and radial distance. 
Based on the results of Lukasiak et al. (1993) the following formula is used to describe 
the dependence of the isotopic ratio (near-Earth) on the solar cycle including a 
dependence on the modulation parameter <£(MV) 
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(3) f(A,Z) 

f Source (A,Z) + [exp(yVd>) - 1] 


where y is found by interpolating from the source and near solar maximum values. The 
second term in Eq.(3) can be positive or negative as dependent on the conditions for a net 
gain or reduction in the cosmic nuclei during transit from the galactic sources (Fields et 
al., 1994) to the inner heliosphere. In Table-1 we have used data on isotopic fractions 
near solar maximum where the isotopic fraction for nuclei produced within the 
heliosphere are expected to be at a maximum (Hesse et al., 1991; Lukasiak et al., 1993; 
Webber et al., 1985, 1990; Wiedenback et al, 1981, 1985). For several elements listed in 
Table-1, information on solar modulation was not available and isotopic fractions for 
near-Earth and GCR source were set identical using the near-Earth estimate. 

We represent the primary (near-Earth) 2 H and 3 He spectra as 

(4a) Ke( E ) = {0.0764 + 0.097 exp[-0.5(ln(E / 1 660) / 1 .306) 2 (E) 

(4b) 0 2 h (E) =0.2 <f> iHe (E) 

where foHe(E) and (]) 2 h(E) are subtracted from the Badhwar and O’Neill model (1992) for 
(j)4He(E) and <j),H (E), respectively. Figure-1 shows the assumed dependence of the isotope 
abundance scaling as a function on the solar modulation parameter described by equation 

1 o i 

(3) for C and N. Similar results are found for other primary isotopes listed in Table-2. 
Examples of the GCR energy spectra for hydrogen and helium are shown in Figure-2a 
and for the isotopes of Ne, Si, and Fe at solar minimum in Figure-2b. 

ISOTOPIC EFFECTS IN GCR TRANSPORT 

The reduction of the full 3-dimensional Boltzmann transport equation within the straight- 
ahead and continuous slowing down approximations and assuming velocity conservative 
fragmentation events accurately describes the transport of GCR heavy ions (Wilson et al., 
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1986, 1991, 1995). In this model the heavy ion flux, </)j(E,x) of an ion j with mass number 
Aj, charge number Z A energy E (in units of MeV/u) at shielding depth x (in units of g/cm ) 
is determined by the partial differential equation (Wilson et al., 1991) 


(5) /-/—/ + a^EW^x.E) = (£)&(£,*) 

ox Sj ( E ) dE k 

where a/E) is the energy-dependent absorption cross-section (cm' 1 ) and csy-(E) is the 
fragmentation cross-section for producing an ion j from k. The solution to eq. (5) in the 
HZETRN code is found using the methods of characteristics where the coordinate 
transformation 


(6a) Jjj = x - Rj (E ) ; = x + R/E) 

and the scaled flux 

(6b) ZjOlj'tj) = (E)0j (x, E) 


are introduced leading to the transport equation 


(V) 



+ cr. 




Yj^j-k—xink^k) 


where S/E) is the stopping power of ion j, y, = Zf/Aj, oj the total absorption cross- 
section, and ojk the fragmentation cross-section to produce isotope j from projectile 
isotope k. The numerical solution to Eq.(7) uses a marching procedure (Wilson et al. 
(1990)). The solution of the transport equation for light ions and neutrons is distinct from 
Eq.(7) because of the broad re-distribution of energy in collision events and is described 
elsewhere (Wilson et al., 1991, Cucinotta et al., 1995b, and Clowdsley et al., 2000). 

Because Eq.(7) is a coupled integro-differential equation for the many GCR primary and 
secondary nuclei, required computer memory allocations increase rapidly as the number 
of ion species is increased and was an important consideration in the early 1990’s. 



However, such practical limitations no longer exist, even on small computer workstations 
with sufficient RAM (e.g. Pentium III or higher with 0.5 Giga-bytes or more), and an 
unrestricted isotopic-grid can presently be implemented with no memory or storage 
problems. When using a reduced isotopic grid, fragments not contained in the grid are 
assigned to a nearby mass of identical charge number. This introduces two types of 
errors: First the range-energy and stopping powers are altered from their true values. This 
error is expected to be small at high-energies for A» 1 since here ion ranges are well 
described by A/Zf scaling factors. The second type of error occurs in the absorption and 
fragmentation cross-sections. Here because of the re-assignment of the mass number 
when using a reduced grid, an error is introduced by the change in neutron number from 
its true values. The fluence spectra of high-energy neutrons produced in fragmentation 
events is also modified by forcing the physics onto a reduced isotopic-grid. This latter 
error is expected to be small for heavy target nuclei (^4> 1 6) because neutron production is 
dominated by light-particle ( n , p, d, t, h, and a) interactions on target nuclei, but may be 
non-negligible for light target atoms (A< 16) (Cucinotta et al. 1998a). Since materials with 
high-hydrogen content are known to be the optimal shielding materials, the changes in 
the neutron fluence due to the use of a full isotopic-grid should be considered. 

QUANTUM FRAGMENTATION MODEL 

We next discuss the quantum multiple scattering description of heavy ion fragmentation 
(QMSFRG), which has been quite successful in describing the physics of the abrasion - 
ablation model of fragmentation and experimental data (Cucinotta, et al. 1992, 1994a, 
1997, 1998b). For inclusive reactions where a single fragment originating in the 
projectile is measured, closure is performed on the final target state with a momentum 
vector denoted p v used to represent these states. The total momentum transfer is 
q = Pr - Px where p T is the initial target momentum. The pre-fragment, F* excitation 
spectrum following nucleon or alpha particle abrasion can be represented in terms of an 
impact parameter dependent convolution of the pre-fragment excitation response for a 
transition of the pre- fragment core from state h to «’ and the project fireball response 
(Cucinotta et al., 1994a, 1998b) 
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( 8 ) 


da 
dE F * 


< T | j d 2 qd 2 bd 2 b' e !q(b - b , p^ (6, b')A n n , (q, b, b, E F , ) | T > 


where b (b ) is the impact parameter, and q the momentum transfer. The abrasion 
response is defined as the interaction of the projectile fireball with the target after 
performing closure over the final fireball states 


( 9 ) 


A nA^b’Af) 


r die 

j — Sj. < R'\ Q„\b') \RxR\ Q„{b) I R' > S(E, 

1 (2x) 



where the Qrt represent the fireball-target profile operator, and kf< the projectile fireball 
momentum vector. The abrasion-response represents a complicated many-body operator 
that is solved by approximation using closure over the target and fireball states for 
evaluating the pre-fragment distribution. The one-particle abrasion response has been 
evaluated using the shell model response functions (Cucinotta and Dubey, 1993). The 
pre- fragment excitation is described in terms of the transition matrix 

(10) P ny (b,b',E F ,) = <F\\ Q F , T \b') | F *„, >< F *„, | Q F * T (b) | F\> 


where matrix elements for the pre-fragment excitation are evaluated over the many-body 
profile operators, Qft ■ In the model a convolution approach is used to derive the mutli- 
knockout spectrum from the single-fragmentation term (Cucinotta and Dubey, 1994a). 
The resulting excitation spectrum is broad with a shape similar to a log-normal 
distribution with mean energies from 20 to 30 MeV for one-nucleon removal (Cucinotta 
and Dubey, 1993; 1994a). 

The de-excitation of the pre-fragments in nuclear ablation is described in a stochastic 
process using a Master equation for nuclear de-excitation by particle emission (Cucinotta 
and Wilson, 1996b). If /(E,t) is the probability of finding the nuclei b at time t with 
excitation energy £), and Pk(E) be the probability that the nuclei, b will emit ion k with 
energy E, then the Master equation is 
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(11) dfh{F ;; J) = Xf dEf a {Elt)P;{E)-Y J [dEf b {El,t)P b k {E) 

at j k 

In Eq.(ll) the first-term on the right corresponds to gains by decays a -> b+j and the 
second term from losses due to decays b -> c+k where the j (or k) are light-particle 
emissions ( n , p, d, t, h, or a). Equation (11) is solved by iteration up to excitation 
energies of 150 MeV and by approximation for higher values (Cucinotta et al. 1996b, 
1997a). Important features of this solution is the correct description of the nuclear level- 
density including nuclear-shell effects at low-excitation energies, and the use of measured 
values for the nuclear masses. The fragmentation cross-section is then evaluated from as 

(12) cr F (A F ,Z F ) = £f dE F > ^ f(A F „Z F . >A f ,Z f ) 

where / (A F , , Z F , > A F ,Z F ) is the solution to Eq.(l 1). 

The iso-spin dependence in fragmentation cross-sections enters in several ways. First, in 
the many-body profile functions, the energy-dependent two-body nucleon interaction 
parameters are summed over the possible projectile and target nucleon scattering 
combinations. Second, the nuclear wave-functions are dependent on shell structure. The 
largest contribution to iso-spin effects occur in the nuclear ablation process (de- 
excitation) where the pre-fragments formed and their level spectra are greatly influenced 
by the projectile or pre-fragment iso-spin. The model uses Coulomb trajectories 
(Cucinotta et al., 1997a) and in-medium nucleon-nucleon interaction (Tripathi et al., 
2000), which are important at lower energies (<200 MeV/u). 

RESULTS 

We first illustrate the accuracy of the QMSFRG model and the effects of iso-spin on 
fragmentation cross-sections. Figures 3 and 4 show comparisons of the model to 
experimental data for the elemental distributions of fragments for several nuclei of 
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similar mass number where nearby isotopes are chosen for comparison. The data of 
Webber et al. (1998) used fragmentation of a primary beam and steering magnets to form 
several less abundant nuclei, and we note the beams are only about 90% pure for 29 Si, 
46 Sc, and 46 Ti results shown. In all comparisons the one-nucleon removal cross-sections 
by electro-magnetic dissociation are included (Norbury et al., 1988). The iso-spin, T z =0 
nuclei display large odd-even effects, which are reduced for the T z + 0 nuclei. The odd- 
even effects are present for all target nuclei, however smaller for hydrogen targets due to 
the reduced abrasion probability for large mass removal on hydrogen. One of errors that 
results from transporting ions using a reduced mass-grid is seen by comparing 
fragmentation cross-sections for nearby projectiles where large differences in many of the 
production cross-sections occur for neighboring projectile nuclei. The model accurately 
reproduces the effects observed in the experiments. Figures-5 show comparisons of 
QMSFRG predictions to experimental data for the isotopic distribution of fragments from 
40 Ar and 56 Fe, and Figures 6-8 correlations between the model and experimental data for 
20 Ne, 32 S, 36 Ar, and 40 Ar and 56 Fe elemental production cross sections . The results show 
good agreement between theory and experiments with the QMSFRG model within 25% 
percent of the experimental values for over 90% of the data in the literature. 

For our extension of the HZETRN code, Table-2 shows the isotopic table of 170 nuclei 
developed as a complete list of GCR primary and secondary nuclei to be used in GCR 
transport problems. We have also listed in the Table-2 the iso-spin, and half-life for the 
unstable nuclei along with the decay mode. This table of nuclei includes all nuclei of 
significant abundance with iso-spin +3/2 to -3 that appear in GCR transport problems. By 
contrast early versions of the HZETRN code used a much smaller ion table. The 
expanded grid used here will allow for improved description of the physics, and to 
discuss many applications where unstable nuclei are central to understanding. Also listed 
in the first column of Table-2 is the index scheme that is used in the HZETRN code, 
which is used to label the shielding-depth and energy dependent fluence matrix. The 
index for the nuclei are ordered by increasing mass number, A followed by charge 
number, Z for a given A. 
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The GCR nuclei are completely stripped and therefore the decay mode and half-life for 
unstable nuclei could be differential from those observed in laboratories on Earth. Garcia- 
Munoz et al. (1987) has noted the following differences: 1) electron capture branches, 
which are inactive for GCR nuclei relative to (1-decay, 2) for non-stripped nuclei there 
will be two S-shell electrons that participate in electron capture that will not be available 
in the decay of GCR nuclei, and 3) changes in screening effects. They have made 
estimates of the elongation of the half-time for fully stripped nuclei due to these 
processes, which indicate an approximate doubling of the decay time observed for 
laboratory nuclei that normally decay by electron capture. Since these decay times are 
much longer than the transit time of nuclei in shielding they are not considered here. 
However, it will useful in the future to further consider these processes when studying the 
effects of stopping GCR nuclei on planetary atmospheres or surfaces and in tissues. 

Figure-9 shows results from the HZETRN code at solar minimum ( 0=428 MV) behind 5 
g/cm 2 of aluminum shielding. Comparison of the mass-flux spectra for a 59-isotope grid 
and the 170-isotope grid are shown. The 170-isotope grid was developed by considering 
the fragmentation cross-sections for a large number of GCR primary nuclei and dominant 
fragments in several materials. In Figure-9 we have scaled the flux by the square of the 
ion charge as a measure of the ionization power of each mass group. Large differences 
are seen for many nuclei. Figure-lOa shows the percent error in the mass-fluence spectra 
resulting from the use of the reduced grid for depths of 5 and 20 g/cm of aluminum. 
Errors greater than 100% are seen for many nuclei, however in most cases such large 
errors only occur for the less abundant nuclei. The probability of biological effects are 
expected to increase in a ma nn er proportional to Z 2 for a given energy, and the elemental- 
flux distribution may be a sufficient test of transport models for supporting exploration 
studies. In Figure-lOb we show a similar comparison to that of Figure-lOb, however 
here for the elemental-fluence spectra. The errors are indeed less substantial than those of 
the mass-fluence spectra, yet are larger than 10% in several cases. Similar comparisons 
near solar maximum conditions are shown in Figures-11. In Table-3 results for the 
elemental and neutron excess dependence of the point dose equivalent behind aluminum 
shielding are shown. The 7<-l nuclei are not significant, while all other cases make 
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important contributions. Table-4 shows the fluence at several depths for several of the 
cosmic-ray clock nuclei and other less abundant nuclei of interest for scientific studies. 


DISCUSSION 

It has been recognized for many years that for the description of GCR transport in 
shielding, theoretical models and experimental data describing the nuclear interactions 
and propagation of protons, heavy ions and their secondaries, leading to accurate and 
computational efficient transport codes are needed. In the last 25 years such descriptions 
have improved dramatically with major milestones that include the development of an 
accurate free space GCR model (Badhwar and O’Neill, 1992), the HZETRN code 
(Wilson, 1977; Wilson, et al., 1990), the measurement of a significant number of 
fragmentation cross-sections (Brechtmann and Heinrech, 1988; Webber et al., 1990b 
Knott et al. 1996, 1997; and Zeitlin et al., 1997, 2001, 2002), and the development of an 
accurate nuclear fragmentation model (Cucinotta et al., 1998a). Laboratory (Shimmerling 
et al, 1989), and spaceflight (Badhwar and Cucinotta, 2000) validation data have also 
become available during this time period. The combination of the GCR model of 
Badhwar and O’Neill, QMSFRG cross section data base, and HZETRN transport code 
have been shown to agree with flight measurements of GCR dose and dose equivalent 
within +15% on several space vehicles. However, further spectral data sets, both in space 
and at heavy ion accelerators, are needed to fully validate these codes. 

The implementation of heavy ion transport models has progressed from models that did 
not satisfy unitarity (Letaw et al., 1983), to the current fully energy-dependent models 
with accurate absorption cross-sections (Shinn et al., 1994; Wilson et al., 1993; 
Cucinotta, 1993). Future work may still be required for light-particle transport ( n , p, d, t, 
h, a, and mesons and their decays), including establishing production cross-section 
models and data, and understanding the role of angular deflections, which are more 
important for neutrons than for heavy ion transport. However, the heavy ion problem is in 
much better shape with many of the remaining task ones of implementation. One 
exception may be improvements in fragmentation cross-sections and laboratory 
validation for the Z=1 to 5 nuclei produced from the heavier projectile nuclei (Z>10). 
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The present paper addressed two implementation tasks; the use of a free-space GCR 
model, which includes the isotopic composition of the primaries, and the extension of the 
HZETRN code to a complete isotopic-grid. The problem of the isotopic distribution of 
the primary GCR and their modulation during the solar cycle has been treated in a 
parametric way in this paper. In this approach, we maintain the accuracy of the GCR 
modulation model for the GCR elemental spectra from Badhwar and O’Neill (1992), 
however re-distribute the fluence of each element amongst its isotopes using estimates 
from satellite data. The coupling of the HZETRN code to the Leaky Box Model is a 
possible approach to a more exact treatment of this problem including the description of 
the energy dependence of the isotopic fractions. However, for shielding applications 
these effects are expected to be small based on the current study. The use of the complete 
isotopic grid of nuclei allows HZETRN applications on the study of the so-called cosmic 
ray “clocks” with lifetimes similar to the time spent by GCR nuclei in the galaxy (~1 M- 
yr) as well as studies of signature nuclei from the decay of GCR nuclei in the Mars 
atmosphere or planetary surfaces. Because the use of a reduced-grid leads to error and 
there are no practical limitations in using a complete isotopic-grid at this time, we 
recommend it be used when initiatives to design space exploration vehicles begin. Future 
tasks that remain are to implement physical models of the GCR isotopic environment and 
to continue to refine the QMSFRG model including comparisons to new fragmentation 
data as they become available. For the many deformed and highly deformed projectile 
nuclei considered herein, methods to consider this deformation in the abrasion process 
are needed. Other considerations for future work is to study the isotopic effects for other 
materials and to consider the decay of radioactive isotopes produced as GCR secondaries 
in the Mars atmosphere, shielding materials or tissue. 
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Table la: Isotopic Composition of GCR Elements Z= 3 to 12. 


Isotope 

Near-Earth Fraction 

Source Fraction 


Z=3 


6 Li* 

0.5 

0.5 

7 Li 

0.5 

0.5 


Z=4 


7 Be* 

0.5 

0.5 

9 Be 

0.35 

0.35 

10 Be 

0.15 

0.15 


Z=5 


10 B 

0.31 

0.2 

n B 

0.69 

0.8 


Z=6 


12 C 

0.92 

0.999 

13 c 

0.08 

0.001 


Z=7 


14 n 

0.48 

0.78 

15 n 

0.57 

0.22 


Z=8 


16 0 

0.946 

0.985 

17 0 

0.027 

0.008 

18 0 

0.027 

0.007 


Z=70 


20 Ne 

0.55 

0.68 

21 Ne 

0.10 

0.0 

22 Ne 

0.35 

0.32 


Z= 72 


24 Mg 

0.64 

0.74 

25 Mg 

0.18 

0.14 

26 Mg 

0.18 

0.13 


*Data on solar modulation was not found and thus near-Earth and source composition are set 
equal. 






Table lb: Isotopic Composition of GCR Elements Z= 13 to 20. 


Isotope 

Near-Earth Fraction 

Source Fraction 


Z=13 


26 A1 

0.02 

0.0 

27 ai 

0.98 

1.0 


Z = 14 


28 Si 

0.84 

0.902 

29 Si 

0.08 

0.054 

30 Si 

0.08 

0.044 


Z=16 


32 S 

0.69 

0.96 

33 s 

0.15 

0.02 

34 s 

0.16 

0.02 


Z=17 


35 C1 

0.52 

1.0 

36 C1 

0.41 

0.0 

37 C1 

0.26 

0.0 


Z=18 


36 Ar 

0.64 

1.0 

37 At 

0.03 

0.0 

38 Ar 

0.30 

0.0 

5-h 

< 

o 

"St 

0.03 

0.0 


Z=20 


40 Ca 

0.4 

1.0 

41 Ca 

0.2 

0.0 

42 Ca 

0.2 

0.0 

43 Ca 

0.2 

0.0 

44 Ca 

0.2 

0.0 
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Table lc: Isotopic Composition of GCR Elements Z=22 to 26. 


Isotope 

Near-Earth Fraction 

Source Fraction 


Z=22 


46r jq* 

0.27 

0.27 

47 Ti 

0.31 

0.31 

48 Ti 

0.30 

0.30 

49 Ti 

0.09 

0.09 

5° T i 

0.03 

0.03 


Z=23 


49 v* 

0.53 

0.53 

50y 

0.24 

0.24 

51 v 

0.23 

0.23 


Z=24 


50 Cr* 

0.26 

0.26 

51 Cr 

0.26 

0.26 

52 Cr 

0.48 

0.48 


Z=25 


53 Mn* 

0.43 

0.43 

54 Mn 

0.17 

0.17 

55 Mn 

0.40 

0.40 


Z=26 


54 Fe 

0.076 

0.055 

55 Fe 

0.084 

0.078 

56 Fe 

0.763 

0.792 

57 Fe 

0.076 

0.075 


*Data on solar modulation was not found and thus near-Earth and source composition are set 
equal. 
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Table 2a. Isotopic Grid Used in Present HZETRN Calculations for Elements Z=0-1 1. 


Index, j 

Nuclei 

Tz 

Lifetime 

Decay 

Index, j 

Nuclei 

Tz 

Lifetime 

Decay 

1 

'n 

-1/2 

614.6 s 

p- 

24 

T2 n 

+ 1 

0.0111 s 

e+, p+ 

2 

’H 

+1/2 

Stable 


27 

13 n 

+1/2 

9.965 m 

e+, p+ 

3 

2 H 

0 

Stable 


29 

14 n 

0 

Stable 


4 

3 h 

-1/2 

12.33 y 

P- 

32 

15 n 

- 1/2 

Stable 


5 

3 He 

+1/2 

Stable 


34 

16 n 

-1 

7.13 s 

P- 

6 

4 He 

0 

Stable 


36 

17 n 

-3/2 

4.173 s 

P- 

7 

6 He 

-1 

0.807 s 

P- 

30 

14 0 

+ 1 

70.606 s 

e+, p+ 

8 

6 Li 

0 

Stable 


33 

15 o 

+1/2 

122.24 s 

e+, p+ 

9 

7 Li 

-1/2 

Stable 


35 

16 o 

0 

Stable 


11 

8 Li 

-1 

0.838 s 

P- 

37 

17 0 

- 1/2 

Stable 


13 

9 Li 

-3/2 

0.178 s 

P- 

39 

18 0 

-1 

Stable 


10 

7 Be 

+1/2 

53.12 d 

e 

42 

19 0 

- 3/2 

26.91 s 

P- 

14 

9 Be 

-1/2 

Stable 


45 

20 o 

-2 

13.51 s 

P- 

16 

10 Be 

-1 

1.51x10® y 

p- 

38 

17 f 

+1/2 

64.49 s 

e+, p+ 

19 

^Be 

- 3/2 

13.81 s 

p- 

40 

18 f 

0 

0.1098 s 

e+, p+ 

12 

8 B 

+1 

0.770 s 

p+ 

43 

19 f 

- 1/2 

Stable 


15 

9 b 

+1/2 

0.54 keV 

2ap 

46 

20 F 

-1 

11.0s 

P- 

17 

10 B 

0 

Stable 


48 

21 F 

-3/2 

4.158 s 

P- 

20 

"B 

-1/2 

Stable 


41 

18 Ne 

+1 

1 .672 s 

e+, p+ 

22 

12 b 

-1 

0.0202 s 

P- 

44 

19 Ne 

+1/2 

17.22 s 

e+, p+ 

25 

13 b 

- 3/2 

0.01736 s 

P- 

47 

20 Ne 

0 

Stable 


18 

10 c 

+1 

19.255 s 

e+, p+ 

49 

21 Ne 

- 1/2 

Stable 


21 

"c 

+1/2 

0.0204 s 

e+, p+ 

51 

22 Ne 

-1 

Stable 


23 

12 c 

0 

Stable 


53 

23 Ne 

- 3/2 

37.24 s 

P- 

26 

13 c 

-1/2 

Stable 


56 

24 Ne 

-2 

3.38 m 

P- 

28 

14 c 

-1 

5730 y 

P- 

50 

21 Na 

+1/2 

22.49 s 

e+, p+ 

31 

15 c 

-3/2 

2.449 s 

P- 

52 

22 Na 

0 

2.6019 y 

e+, p+ 






54 

23 Na 

-1/2 

stable 







57 

24 Na 

-1 

14.959 h 

P- 






59 

25 Na 

- 3/2 

59.1 s 

P- 






62 

26 Na 

-2 

1 .072 s 

P- 




Table 2b. Isotopic Grid Used in Present HZETRN Calculations for Elements Z=12-18. 


Index, j 

Nuclei 

Tz 

Lifetime 

Decay 

Index, j 

Nuclei 

Tz 

Lifetime 

Decay 

55 

ZJ Mg 

+1/2 

11.317s 

e+, p+ 

78 

30g 

+1 

1.178s 

e+, p+ 

58 

24 Mg 

0 

stable 


81 

31 s 

+1/2 

2.572 s 

e+, p+ 

60 

“Mg 

-1/2 

stable 


84 

32 s 

0 

Stable 


63 

Mg 

-1 

stable 


87 

33 s 

- 1/2 

Stable 


66 

27 Mg 

-3/2 

9.458 m 

p- 

90 


-1 

Stable 


69 

ho 

00 

CQ 

-2 

20.91 h 

p- 

93 

35 s 

-3/2 

87.32 d 

P- 

61 

25 ai 

+1/2 

7.183 s 

e+, p+ 

96 

36 s 

-2 

Stable 


64 

26 ai 

0 

7.17x10 s y 

e+, p+ 

99 

37 s 

- 5/2 

5.05 m 

P- 

67 

27 ai 

-1/2 

stable 


103 

38 s 

-2 

170.3 m 

P- 

70 

28 ai 

-1 

2.241 m 

P- 

88 

33 CI 

+1/2 

2.511 s 

e+, p+ 

72 

29 ai 

- 3/2 

6.56 m 

P- 

91 

34 CI 

0 

1 .5264 s 

e+, p+ 

75 

30 ai 

-2 

3.60 s 

P- 

94 

35 CI 

- 1/2 

Stable 


65 

26 Si 

+1 

2.234 s 

e+, p+ 

97 

36 CI 

-1 

3.01x10 s y 

P- 

68 

27 Si 

+1/2 

4.16s 

e+, p+ 

100 

37 CI 

-3/2 

Stable 


71 

28 Si 

0 

stable 


104 

38 CI 

-2 

37.34 m 

P- 

73 

29 Si 

-1/2 

stable 


108 

39 CI 

- 5/2 

55.60 m 

P- 

76 

30 Si 

-1 

stable 


95 

35 Ar 

+1/2 

1.775 s 

e+, p+ 

79 

31 Si 

- 3/2 

0.1573 s 

P- 

98 

36 Ar 

0 

Stable 


82 

32 Si 

-2 

150 y 

P- 

101 

37 Ar 

- 1/2 

35.04 d 

e 

85 

33 Si 

- 5/2 

6.18s 

P- 

105 

38 Ar 

-1 

Stable 


74 

29 p 

+1/2 

4.140 s 

e+, p+ 

109 

39 Ar 

- 3/2 

269 y 

P- 

77 

30 p 

0 

2.498 m 

e+, p+ 

112 

40 Ar 

-2 

Stable 


80 

31 p 

-1/2 

stable 


115 

41 Ar 

-5/2 

109.34 m 

P- 

83 

32 p 

-1 

14.262 d 

P- 

118 

42 Ar 

-3 

32.9 y 

P- 

86 

33 p 

- 3/2 

25.34 d 

P- 






89 

34 p 

-2 

12.43 s 

P- 






92 

35 p 

- 5/2 

47.3 s 

P- 
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Table 2c. Isotopic Grid Used in Present HZETRN Calculations for Elements Z=19-28. 


Index, j 

Nuclei 

Tz 

Lifetime 

Decay 

Index, j 

Nuclei 

Tz 

Lifetime 

Decay 

102 

37 R 

+1/2 

1 .226 s 

e+, p+ 

134 


0 

0.422 s 

e+, p+ 

106 

38 k 

0 

7.636 m 

e+, p+ 

137 

47 v 

- 1/2 

32.6 m 

e+, p+ 

110 

39 k 

- 1/2 

Stable 


140 

48 v 

-1 

15.9735 d 

e+, p+ 

113 

40 k 

-1 

1. 277x1 0 9 y 

P- 

143 

49 v 

- 3/2 

330 d 

e 

116 

41 K 

- 3/2 

Stable 


146 

5°V 

-2 

1.4x10 17 y 

CD 

+ 

CO 

+ 

119 

42 k 

-2 

12.36 h 

P- 

149 

51 v 

- 5/2 

Stable 


122 

43 k 

- 5/2 

22.3 h 

P- 

152 

52 v 

-3 

3.743 m 

p- 

107 

38 Ca 

+ 1 

0.440 s 

e+, p+ 

141 

48 Cr 

0 

21.56 h 

e+, p+ 

111 

39 Ca 

+1/2 

0.8596 s 

e+, p+ 

144 

49 Cr 

- 1/2 

42.3 m 

e+, p+ 

114 

40 Ca 

0 

Stable 


147 

50 Cr 

-1 

>1 .8x10 17 y 

e 

117 

41 Ca 

-1/2 

1.03x10 5 y 

e 

150 

51 Cr 

-3/2 

27.7025 d 

e 

120 

42 Ca 

-1 

Stable 


153 

52 Cr 

-2 

Stable 


123 

43 Ca 

- 3/2 

Stable 


155 

53 Cr 

- 5/2 

Stable 


125 

44 Ca 

-2 

Stable 


158 

54 Cr 

-3 

Stable 


128 

45 Ca 

- 5/2 

162.61 d 

p- 

148 

50 Mn 

0 

0.284 s 

e+, p+ 

131 

46 Ca 

-3 

Stable 


151 

51 Mn 

- 1/2 

0.0462 s 

e+, p+ 

121 

42 Sc 

0 

0.681 s 

e+, p+ 

154 

52 Mn 

-1 

5.591 d 

e+, p+ 

124 

43 Sc 

- 1/2 

3.891 h 

e+, p+ 

156 

53 Mn 

-3/2 

3.74x10® 

e 

126 

44 Sc 

-1 

3.927 h 

e+,p+ 

159 

M Mn 

-2 

312.3 d 

0 

+ 

■co 

+ 

129 

45 Sc 

- 3/2 

stable 


161 

55 Mn 

-5/2 

Stable 


132 

46 Sc 

-2 

83.79 d 

P- 

157 

53 Fe 

- 1/2 

0.00851 s 

0 

+ 

co 

+ 

135 

47 Sc 

- 5/2 

3.3492 d 

P- 

160 

54 Fe 

-1 

Stable 


138 

48 Sc 

-3 

43.67 h 

P- 

162 

55 Fe 

-3/2 

2.73 y 

e 

127 

44 Ti 

0 

63 y 

e 

164 

56 Fe 

-2 

Stable 


130 

45 Ti 

- 1/2 

184.8 m 

e+, p+ 

167 

57 Fe 

- 5/2 

Stable 


133 

46 Ti 

-1 

Stable 


163 

55 Co 

- 1/2 

17.53 h 

e+,p+ 

136 

47 Ti 

- 3/2 

Stable 


165 

56 Co 

-1 

77.27 d 

e+, p+ 

139 

48 Ti 

-2 

Stable 


168 

57 Co 

-3/2 

271.79 d 

e 

142 

49 Ti 

- 5/2 

Stable 


166 

56 Ni 

0 

6.077 d 

e+, p+ 

145 

SOjj 

-3 

Stable 


169 

57 Ni 

- 1/2 

35.60 h 

e+, p+ 






170 

58 Ni 

-1 

Stable 
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Table 3a: Elemental (Z) and Neutron excess (Y) dependence on GCR dose equivalent 
behind 0 g/cm 2 of Aluminum Shielding. 


GCR Dose Equivalent per Year for near Solar Minimum 


z 

Y<0 

Y=0 

Y=1 

Y=2 

Y=3 

Y>3 

Total-Z 

0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1 

9.15 

0.21 

0.00 

0.00 

0.00 

0.00 

9.37 

2 

0.50 

3.61 

0.00 

0.00 

0.00 

0.00 

4.11 

3 

0.00 

0.04 

0.04 

0.00 

0.00 

0.00 

0.08 

4 

0.05 

0.00 

0.03 

0.01 

0.00 

0.00 

0.10 

5 

0.00 

0.23 

0.36 

0.00 

0.00 

0.00 

0.60 

6 

0.00 

3.72 

0.31 

0.00 

0.00 

0.00 

4.03 

7 

0.00 

0.74 

0.99 

0.00 

0.00 

0.00 

1.73 

8 

0.00 

11.26 

0.32 

0.32 

0.00 

0.00 

11.90 

9 

0.00 

0.00 

0.38 

0.00 

0.00 

0.00 

0.38 

10 

0.00 

2.33 

0.45 

1.48 

0.00 

0.00 

4.26 

11 

0.00 

0.00 

1.33 

0.00 

0.00 

0.00 

1.33 

12 

0.00 

6.89 

1.86 

1.88 

0.00 

0.00 

10.63 

13 

0.00 

0.12 

2.34 

0.00 

0.00 

0.00 

2.46 

14 

0.00 

11.22 

1.08 

1.06 

0.00 

0.00 

13.36 

15 

0.00 

0.00 

0.85 

0.00 

0.00 

0.00 

0.85 

16 

0.00 

2.87 

0.59 

0.65 

0.00 

0.00 

4.11 

17 

0.00 

0.00 

0.52 

0.41 

0.26 

0.00 

1.19 

18 

0.00 

1.37 

0.06 

0.68 

0.00 

0.05 

2.15 

19 

0.00 

0.00 

1.71 

0.00 

0.00 

0.00 

1.71 

20 

0.00 

1.57 

0.00 

0.78 

0.78 

0.78 

3.92 

21 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 

1.00 

22 

0.00 

0.00 

0.00 

1.02 

1.17 

1.59 

3.78 

23 

0.00 

0.00 

0.00 

0.00 

1.00 

0.89 

1.88 

24 

0.00 

0.00 

0.00 

1.02 

1.02 

1.89 

3.93 

25 

0.00 

0.00 

0.00 

0.00 

1.26 

1.67 

2.92 

26 

0.00 

0.00 

0.00 

2.05 

2.26 

22.50 

26.80 

27 

0.00 

0.00 

0.00 

0.00 

0.10 

0.00 

0.10 

28 

0.00 

0.00 

0.00 

1.45 

0.00 

0.00 

1.45 

Total-Y 

9.70 

46.18 

13.22 

12.81 

8.85 

29.35 

120.11 







Table 3b: Elemental (Z) and Neutron excess (Y) dependence on GCR dose equivalent 
behind 5 g/cm 2 of Aluminum Shielding. 


GCR Dose Equivalent per Year near Solar Minimum 

Z 

Y<0 

Y=0 

Y=1 

Y=2 

Y=3 

Y>3 

Total-Z 

0 

0.00 

0.00 

0.00 

0.65 

0.00 

0.00 

0.65 

1 

15.03 

2.56 

1.74 

0.00 

0.00 

0.00 

19.33 

2 

2.48 

15.60 

0.00 

0.00 

0.00 

0.00 

18.08 

3 

0.00 

0.04 

0.04 

0.00 

0.00 

0.00 

0.09 

4 

0.05 

0.00 

0.04 

0.02 

0.00 

0.00 

0.12 

5 

0.00 

0.23 

0.34 

0.00 

0.00 

0.00 

0.58 

6 

0.03 

3.14 

0.29 

0.01 

0.00 

0.00 

3.47 

7 

0.01 

0.62 

0.83 

0.01 

0.00 

0.00 

1.48 

8 

0.07 

8.67 

0.27 

0.27 

0.00 

0.00 

9.28 

9 

0.01 

0.03 

0.31 

0.01 

0.01 

0.00 

0.37 

10 

0.03 

1.69 

0.37 

1.10 

0.01 

0.00 

3.20 

11 

0.01 

0.06 

0.99 

0.04 

0.01 

0.00 

1.11 

12 

0.05 

4.71 

1.33 

1.34 

0.01 

0.00 

7.44 

13 

0.01 

0.15 

1.66 

0.03 

0.01 

0.00 

1.87 

14 

0.06 

7.46 

0.77 

0.75 

0.01 

0.00 

9.06 

15 

0.00 

0.02 

0.60 

0.04 

0.02 

0.01 

0.70 

16 

0.02 

1.86 

0.42 

0.48 

0.02 

0.01 

2.82 

17 

0.00 

0.02 

0.38 

0.32 

0.21 

0.01 

0.93 

18 

0.01 

0.89 

0.08 

0.51 

0.04 

0.06 

1.59 

19 

0.00 

0.02 

1.12 

0.04 

0.04 

0.04 

1.27 

20 

0.01 

1.03 

0.03 

0.58 

0.57 

0.58 

2.81 

21 

0.00 

0.00 

0.02 

0.05 

0.70 

0.08 

0.85 

22 

0.00 

0.01 

0.03 

0.72 

0.81 

1.12 

2.69 

23 

0.00 

0.00 

0.02 

0.05 

0.71 

0.66 

1.44 

24 

0.00 

0.00 

0.03 

0.68 

0.73 

1.33 

2.77 

25 

0.00 

0.00 

0.02 

0.05 

0.88 

1.23 

2.18 

26 

0.00 

0.00 

0.03 

1.34 

1.50 

13.90 

16.77 

27 

0.00 

0.00 

0.00 

0.01 

0.06 

0.00 

0.08 

28 

0.00 

0.00 

0.01 

0.87 

0.00 

0.00 

0.88 

Total-Y 

17.88 

48.84 

11.80 

9.95 

6.37 

19.03 

113.88 
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Table 3c: Elemental (Z) and Neutron excess (Y) dependence on GCR dose equivalent 
behind 20 g/cm 2 of Aluminum Shielding. 


GCR Dose Equivalent per Year near Solar Minimum 

z 

Y<0 

Y=0 

Y=1 

Y=2 

Y=3 

Y>3 

Total-Z 

0 

0.00 

0.00 

0.00 

2.00 

0.00 

0.00 

2.00 

1 

18.95 

3.09 

2.21 

0.00 

0.00 

0.00 

24.24 

2 

2.67 

17.03 

0.00 

0.00 

0.00 

0.00 

19.69 

3 

0.00 

0.03 

0.03 

0.00 

0.00 

0.00 

0.07 

4 

0.04 

0.00 

0.04 

0.02 

0.00 

0.00 

0.10 

5 

0.00 

0.13 

0.21 

0.00 

0.00 

0.00 

0.35 

6 

0.05 

1.51 

0.13 

0.02 

0.00 

0.00 

1.72 

7 

0.02 

0.41 

0.38 

0.01 

0.00 

0.00 

0.83 

8 

0.11 

3.69 

0.12 

0.11 

0.00 

0.00 

4.02 

9 

0.01 

0.04 

0.17 

0.02 

0.01 

0.00 

0.26 

10 

0.04 

0.75 

0.16 

0.46 

0.01 

0.00 

1.43 

11 

0.01 

0.09 

0.45 

0.05 

0.02 

0.00 

0.62 

12 

0.08 

1.84 

0.52 

0.50 

0.01 

0.00 

2.95 

13 

0.02 

0.12 

0.69 

0.04 

0.02 

0.00 

0.89 

14 

0.08 

2.57 

0.29 

0.26 

0.02 

0.01 

3.22 

15 

0.01 

0.04 

0.26 

0.05 

0.03 

0.01 

0.39 

16 

0.03 

0.71 

0.15 

0.18 

0.03 

0.02 

1.12 

17 

0.00 

0.02 

0.19 

0.12 

0.08 

0.02 

0.44 

18 

0.01 

0.34 

0.07 

0.19 

0.05 

0.04 

0.70 

19 

0.00 

0.03 

0.39 

0.06 

0.05 

0.05 

0.58 

20 

0.02 

0.47 

0.05 

0.20 

0.19 

0.19 

1.11 

21 

0.00 

0.01 

0.03 

0.07 

0.29 

0.10 

0.49 

22 

0.00 

0.01 

0.04 

0.30 

0.33 

0.46 

1.14 

23 

0.00 

0.00 

0.03 

0.06 

0.30 

0.31 

0.71 

24 

0.00 

0.01 

0.04 

0.26 

0.31 

0.55 

1.16 

25 

0.00 

0.00 

0.02 

0.06 

0.36 

0.55 

1.00 

26 

0.00 

0.00 

0.03 

0.43 

0.53 

4.13 

5.13 

27 

0.00 

0.00 

0.00 

0.01 

0.02 

0.00 

0.04 

28 

0.00 

0.00 

0.01 

0.24 

0.00 

0.00 

0.25 

Total-Y 

22.15 

32.94 

6.69 

5.73 

2.68 

6.45 

76.64 
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Table 4: Annual fluence of several unstable GCR nuclei versus depth of aluminum 
shielding near solar minimum ( 0=428 MV). 




Fluence 

per cm 2 per year 


Nuclei 

0, g/cm 2 

2, g/cm 2 

5, g/cm 2 

20, g/cm 2 

6 He 

0.0 

9-lxlO 1 

2.1x10 2 

5.6x1 0 2 

8 Li 

0.0 

1 .0x1 0 3 

2.2x1 0 3 

5.4x1 0 3 

9 Li 

0.0 

3.0x1 0 1 

7.2x1 0 1 

1 .8x1 0 2 

10 Be 

8.8x1 0 3 

9.1x10 3 

9.3x1 0 3 

8.8x1 0 3 

14 C 

0.0 

7.1x10 2 

1 .7x1 0 3 

3.4x1 0 3 

18 o 

9.9x10 3 

9.3x1 0 3 

8.4x1 0 3 

5.1x10 3 

26 ai 

5.7x1 0 2 

7.6x1 0 2 

9.4x1 0 2 

1 .0x1 0 3 

36 CI 

8.1x10 2 

7.7x1 0 2 

7.1x10 3 

4.4x1 0 2 
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Figure 1 : Parametric model for describing the change in isotopic composition with the 

solar modulation parameter, <P(MV). 
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E, MeV/u 


Figure 2a: Energy spectra for hydrogen and helium isotopes near solar minimum 
(cP=428 MV) and solar maximum (<P=1050 MV). 



E, MeV/u 


Figure 2b: Energy spectra for Ne, Si, and Fe isotopes near solar minimum (cZ>=428) 
showing contributions from different isotopes to primary GCR composition. 
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125 


28 


Si (0.508 GeV/u)+ H ~> Z 




Figure 3: Comparisons of the QMSFRG model to experiment for elemental fragment 
distribution for 28 Si and 29 Si on H at 0.51 GeV/u. Experimental data from Webber et al. 
(1998). 
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Figure 4a: Comparisons of the QMSFRG model to experiment for fragmentation of (left) 
40 Ar on H at 0.35 GeV/u and for (right) 40 Ca on H at 0.76 GeV/u (Knott et al. (1996)). 



Figure 4b: Comparisons of the QMSFRG model to experiment for fragmentation of (left) 
46 Sc on H at 0.57 GeV/u and for (right) 46 Ti on H at 0.57 GeV/u (Webber et al. (1998)). 
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Figure 5: Comparisons of the QMSFRG model to experiment for isotopic distribution of fragments for a) 40 Ar on 12 C at 0.6 GeV/u. 
and b) 56 Fe on 12 C interactions at 0.6 GeV/u. Experimental data from Webber et al. (1990b). 
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Figure 6: a) Comparison of QMSFRG to data of Zeitlin et al. (2001 ) for 20 Ne fragmentation on several targets at 0.6 GeV/u. 
And b) comparison to data of Brechtmann et al. (1997) for 32 S fragmentation on several targets at 0.7 or 1 .2 GeV/u. 
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Figure 7: Comparison of QMSFRG to data of lancu et al. (2005) for 36 Ar or 40 Ar fragmentation on several targets at 0.34 GeV/u. 
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Figure 8: Comparison of QMSFRG model to experimental data for 56 Fe fragmentation on several targets at 0.65 GeV/u (Flesch et al. 
(1999)), 1.05 GeV/u (Zeitlin etal. (1997)), and 1.6 GeV/u (Cummings etal., 1990). 



Figure 9: Comparisons of results from the HZETRN Code for the mass flux distribution 
behind 5 g/cm 2 of aluminum shielding for solar minimum conditions comparing transport 
with the reduced 59-isotope grid to transport with a full 170-isotope grid (proton fluence 
shown for A=1). 
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Figure 10a: Comparisons of the error that results from the HZETRN Code for the mass 
fluence distribution near solar minimum when using a reduced 59-isotope grid compared 
to transport with a 170-isotope grid. 
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Figure 10b: Comparisons of the error that results from the HZETRN Code for the 
elemental fluence distribution near solar minimum when using a reduced 59-isotope grid 
compared to transport with a 170-isotope grid. 
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Figure 11a: Comparisons of the error that results from the HZETRN code for the mass 
fluence distribution near solar maximum when using a reduced 59-isotope grid 
compared to transport with al 70-isotope grid. 
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Figure 1 1 b: Comparisons of the error that results from the HZETRN code for the 
elemental fluence distribution near solar maximum when using a reduced 59-isotope grid 
compared to transport with al 70-isotope grid. 
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